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. Random critical points are generically characterized by multifractal properties. In the field of 

' Anderson localization, Mirlin, Fyodorov, Mildenberger and Evers [Phys. Rev. Lett 97, 046803 

, (2006)] have proposed that the singularity spectrum /(a) of eigenfunctions satisfies the exact sym- 

^ . metry f{2d. — a) — f{a) + d — a at any Anderson transition. In the present paper, we analyse the 

Q ' physical origin of this symmetry in relation with the Gallavotti- Cohen fiuctuation relations of large 

I deviation functions that are well-known in the field of non-equilibrium dynamics: the multifractal 

, spectrum of the disordered model corresponds to the large deviation function of the rescaling ex- 

• ponent 7 = (q — d) along a renormalization trajectory in the effective time t = InL. We conclude 

that the symmetry discovered on the specific example of Anderson transitions should actually be 
satisfied at many other random critical points after an appropriate translation. For many-body 
random phase transitions, where the critical properties are usually analyzed in terms of the multi- 
fractal spectrum H{a) and of the moments exponents X{N) of two-point correlation function [A. 
C/3 i Ludwig, Nucl. Phys. B330, 639 (1990)], the symmetry becomes H{2X{1) - a) = H{a) + a- X{1), 

or equivalently A{N) = A(l - iV) for the anomalous parts A{N) = X{N) - NX(1). We present 
numerical tests in favor of this symmetry for the 2D random Q— state Potts model with various Q. 
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INTRODUCTION 



^ , Among the various areas where multifractaUty appears (see for instance [U, B 0, H, H, 0] and references therein), 
' the case of critical points in the presence of frozen disorder is of particular interest. The idea that multifractality 



occurs at criticality has been first established for quantum Anderson localization transitions [8|, |9| and has been 
the subject of very detailed studies (see the reviews [IH)- However multifractaUty is not limited to this type 
of one-particle quantum transitions, but is expected to hold at many random critical points whenever disorder is 
relevant, in particular in disordered s pin models like random ferromagnets [l^, [IB Eli [13, [Hj [3j spin-glasses 
QQ ' or random field spin systems pol . [21I . , as well as in disordered polymer models like directed polymer in random 

■ media [1^ or disordered wetting models [l^]. The only exceptions to these multifractal behaviors seem to be the 

■ "multiscaling" behaviors [1^, which are even stronger than multifractality, that have been found for some critical 
! correlation functions in disordered quantum spin-chains governed by "Infinite disorder fixed points" So the 
' presence of multifractality at criticality seems generic whenever disorder is relevant. 

^\ [ For Anderson localization models, Mirlin, Fyodorov, Mildenberger and Evers [13] have proposed that the singularity 
' spectrum /(a) of critical eigenfunctions satisfies some remarkable exact symmetry at any Anderson transition (see 
^ I more details in section [II] below). In this paper, we analyse the physical origin of this symmetry in relation with the 
well-know symmetry relations of large deviation functions in non-equilibrium dynamical models. We conclude that 
this symmetry is not specific to Anderson transitions but should be satisfied at many other random critical points, 
in particular at many-body disordered phase transitions after an appropriate translation that we describe. As an 
example of application, we consider the two-dimensional random Potts model with various Q and present numerical 
tests of the symmetry concerning the multifractal spectrum of the two-point order parameter correlation functions. 

The paper is organized as follows. In section [ITl we recall the meaning of the symmetry relation for Anderson 
localization models. In section IIII[ we analyse the similarity with the well-known fluctuation relations for large 
deviation functions occurring in the field of non-equilibrium dynamics. This is the central section of the present 
paper. In section HVl we translate the symmetry relation for /(a) in terms of the multifractal spectrum H{a) that 
characterizes the statistics of two-point correlation function at many-body random critical points. In section |Vl we 
present numerical tests in favor of this symmetry for the two-dimensional random Q— state Potts model with various 
Q. Section IVTl contains our conclusions. 
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II. SYMMETRY RELATION AT ANDERSON TRANSITIONS 
A. Reminder on the singularity spectrum f{a) of critical eigenfunctions 

For Anderson localization transitions, the multifractal spectrum f{a) of critical eigenfunctions is defined as follows 
(for more details see the reviews [l^, Eli): m a sample of size L"^, the number Mhip) of points r where the weight 
|'0l(?^)P scales as L^" behaves as 

UL{a) oc (1) 

L — »oo 

The inverse participation ratios, which are the most convenient order parameters of the transition, can be then 
rewritten as an integral over a 

Yq{L)= j d''r1i/.L(f)|2'J ~ /da ~ L-^^"^ (2) 

where the exponent r(g) can be obtained via a saddle-point calculation in a to obtain the Legcndrc transform formula 

[13, [iH 

- T{q) = max [/(a) - qa] (3) 

OL 

The usual normalization condition of eigenfunctions 

dr|^L(rOr = l (4) 

implies Yq^i{L) = 1 so that T{q = 1) = 0. The weight |V'l(?^)P has thus for disorder averaged value 

1 



lV'L(r)P = (5) 



whereas the typical exponent atyp governing the typical decay 



ln|'0L(f)|2 oc -atyp\nL (6) 

L — *oo 

is determined by the maximum f{atyp) = d of the spectrum /(a). These scaling behaviors, which concern individual 
eigenstates ip, can be translated for the local density of states 

PL{E,r) - ^<5(i? - i?„)|V^B„(rOp (7) 

n 

as follows: for large L, when the energy levels become dense, the sum of Eq. [71 scales as 

PL{E,f) <^ L'\^'E{r)\' (8) 
and its moments involve the exponents T{q) introduced in Eq. [2] 

[pL{E,r)]i ^^^j^ with A{q) = T{q) - dig - 1) (9) 

B. Symmetry relation of f{a) 

For any Anderson transition in the so-called 'conventional symmetry classes' Mirlin, Fyodorov, Mildenbcrgcr 
and Evers pTj have proposed that the singularity spectrum /(a) of critical eigenfunctions satisfies the remarkable 
exact symmetry 



f{2d-a) = ,f{a) + d-a 



(10) 



3 



that relates the regions a < d and a > d. (see for instance Fig. 1 of Ref. [33] b for a pictorial representation of 
this symmetry). In terms of the exponents T{q) of Eq. t he symmetry of Eq. 1101 with respect to the value as ^ d 
becomes a symmetry with respect to the value qs — 1/2 [23| 

T{q) - t{1 - q) ^ d{2q - 1) (11) 

In terms of the exponents A(g) of Eq. [9l the symmetry takes the simpler form 

A(q)-A(l-q) (12) 

As discussed in p^ . this type of symmetry has been first derived for non- linear sigma- models p8l |. which are a 
priori only valid for weak disorder. Nevertheless, Mirlin, Fyodorov, Mildenberger and Evers have argued in that 
Eq. [To] should remain exact even at strong disorder as a consequence of universality of critical properties. Eq. 1101 has 
been checked on the exact expansion in d + e up to 4-loop order poj ] , and numerically for various types of Anderson 
transitions, in particular in the Power-Law random banded matrices [27| . in the symplectic 2D Anderson transition 
[sot , i'^ quantum Hall transition [sij, and in the 3D Anderson transition [s^l (the only exception being, to our 
knowledge, the spin quantum Hall transition ^3S] that belongs to the Bogoliubov-de Genncs symmetry class C of the 
symmetry classification [llj) This symmetry has even been measured in recent experiments [33 |. 

In conclusion, the symmetry of Eq. [10] seems very well satisfied at most of Anderson transitions where it has been 
studied. However its physical origin has remained unclear. In particular, an important issue is whether this symmetry 
is specific to critical theories of Anderson transitions, or whether it could be satisfied at other random critical points. 
These questions have been the motivations of the present work. 



III. PHYSICAL ORIGIN OF THE SYMMETRY RELATION 

The multifractal formalism can be seen as a theory of large deviations in the parameter In L (see [35| for a recent 
review on large deviations in statistical physics). It is thus interesting to discuss in this section the similarity that 
exists between the Mirlin-Fyodorov-Mildenberger-Evers symmetry of Anderson transitions described in the previous 
section, and the well-known symmetry relations of large deviation functions that occur in the field of non-equilibrium 
dynamics. 



A. Fluctuations relations of large deviations functions in non-equilibrium dynamical models 

It is clearly impossible to summarize here all the recent developments concerning the various 'fluctuation relations' 
that have been established in the field of non-equilibrium dynamics recently, and we refer the interested reader to 
some recent reviews [s^, 113, H^, H^, 113, EH, 113 and references therein. Here we will simply recall some basic definitions 
that will be useful for our present purposes. The observables Yt which are expected to become extensive in time in the 
large time limit, usually satisfy some large deviation principle: the probability to have a given time-averaged value 
Yf/t = j behaves at large t as 

Prob (^j = ^ cx e*°(^^ (13) 

where G{j) < is called the large-deviation function. The typical value jtyp corresponds to the point where it vanishes 
G{jtyp) = 0. It has been found that in many cases, the large deviation function G(j) satisfies some symmetry relation 
of the form 

GU) = G{-j) + Kj (14) 

where the constant K may contain the physical parameters for the problem at hand (like the reservoirs fugacities, 
the temperature, the external applied field, etc). Equivalently, the generating function of the variable Yf behaves for 
large time t as 



where /J. (A) is the Legendre transform of G{i) 



e^^O = / dje*^^^+"^J^) cx e'^^^^ (15) 



/x(A) =max[Aj+G(j)] (16) 
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Its series expansion in A 



/x(A)= lim -L^=Xjtyp + —a + o{X^) (17) 



yields the successive cumulants of Yt 



Jtyp 



t^+oo t 



. ^ lin. i^^l-i^ (18) 

The symmetry of Eq. [Threads for the generating function /i(A) 

/i(A) = ^l{-K - A) (19) 

The symmetry of Eq. [M] or Eq. [19] is usually called a Gallavotti-Cohcn fluctuation relation. The full domain of 
validity of this type of symmetry in the field of non-equilibrium dynamics is not easy to state, since it has been derived 
in various contexts with different assumptions on the dynamics, which can be either stochastic or deterministic but 
sufficiently chaotic (see the various presentations in the reviews [H, 113, [H, [H, HO, |4l|, H^). Let us remind here the 
most important ideas on the simplest case: for a stochastic dynamics defined by some Markov chain with transition 
probabilities k{C — ^ C") between configurations 

Pt+i(C)==^Pt(C')fc(C"-.C) (20) 
c 

the symmetry relation of large deviation has usually for origin some 'generalized detailed balance relation' of the form 

ky{C C) ^ k-y{C' C)e^y (21) 

where y denotes the increase of the dynamical quantity Yt for the jump C C , and (— y) denotes the increase of the 
dynamical quantity Yt for the jump C" —* C. Eq. [211 implies that a given dynamical trajectory Traj = {Co, Ci, ..Ct} 
characterized by the value YtlTraj) = yi + y2 + ■• + yt, and the reversed trajectory {—Traj) = {Ct, Ct_i, ..Co} 
characterized by the opposite value Yt{—Traj) = —Yt{Traj) have probabilities related by the simple relation 

PjTraj) ^ fcyi(Co Ci)fc^,(Ci C2)...fcyt (Cf„i Ct) ^ KYtiTra,) (22) 

P{~Traj) k.y^{Ct~^Ct^i)...k^y,{C2^Ci)k^y,{Ci^Co) ^ ' 

The generating function of Eq. [12] reads 

(e^^') = ^ P(Traj)e^^'(^™^") = P(-rra7)e(^+^)^*(^™^') (23) 

Traj Traj 

Using the one-to-one change of variable Traj' = —Traj and the antisymmetry relation Yt{—Traj) — —Yt{Traj), one 
obtains 



e 



XY,^ _ \- P(Traj')e-(^+^)^*(^™^') = (e-(^+^)^*) (24) 



Traj' 



which corresponds to the symmetry relation of Eq. [19] for the generating function introduced in Eq. [15] 

To derive a symmetry relation for a large deviation function of some dynamical observable Yt, it is sufficient to justify 
some 'generalized detailed balance relation' of the form of Eq. [21] for the elementary transition probabilities. We refer 
to the reviews [s^, [13, [H, [H, HO, SH, for the description of various physical situations where this can be done 
(system in contact with two reservoirs at different temperatures, system in an external applied field, etc). However, 
as explained in [i^, for any stochastic process there is always a quantity Bt for which it works by construction: it is 
the quantity Bt whose increment b(C, C') is defined by Eq. [5T]with K = 1, i.e. it corresponds to the logarithm of 
the ratio of the two transition rates [43| 
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The functional Bt of the trajectory Traj = {Co,Ci, ..Ct} for which the fluctuation relation holds by construction 
then reads [i^ 

B.iTrai ^ iC„,C..^C,}) ^ ± ,n {'j^^^f^) (26) 

771—1 ^ ^ ^ ' 

We stop here this brief reminder on large deviation functions occurring in non-equilibrium dynamics, and refer the 
interested reader to and to the reviews [H, [13, [H, [H, [13, El, Sll for the physical interpretation in terms of 
irreversibility and entropy production. We now analyse the analogy with multifractal spectra at random critical 
points. 



B. Dictionary of the analogy 

To make the formal analogy between Eq. [M] and Eq. [10] complete, it is convenient to replace the singularity 
exponent a by its distance to d which represents the 'homogeneous' value under rescaling 

j = a-d (27) 

Then the negative function 

5(7) = /(a = d + 7)-d (28) 

satisfies (Eq. ^ 

.9(7) =.9(-7)+7 (29) 

This equation is now exactly similar to Eq. [T3]with K = 1. The relation between the generating function /i(A) of Eq. 
[16] and the exponents A{q) of Eq. [9] reads 

^(A) = max [A7 + .9(7)] 

7 

— max [X{a — d) + f{a) — d] = — t(— A) — d\ ~ d 

a 

= -A(-A) (30) 

so that the symmetry relation of Eq. [19] with A' = 1 is equivalent to Eq. [12] 

Since the physical situations of a non-equilibrium dynamical system and of a static disordered system at criticality 
are a priori completely different, this analogy may seem completely formal at first sight. It is however very suggestive: 
the large time t of the dynamical system corresponds to the logarithm of the large linear scale L for the random 
critical system 

t = lnL (31) 
The 'effective dynamics' for the random system thus corresponds to renormalization process in scale InL. 



C. Stability of the multifractal spectrum upon renormalization 



In critical phenomena, it is well known that critical properties arc stable under coarse-graining. This explains their 
universal character (independence with respect to microscopic details) and why renormalization is an appropriate 
framework. Similarly for random critical points, the multifractal spectrum is expected to be stable under coarse- 
graining, i.e. the formulations we have written above in terms of the microscopic observables (like the eigenfunction 
weight ijp'ir) in Anderson localization models) can be reformulated in terms of coarse-grained observables as follows 
[lot : if the system of volume L"^ is decomposed into boxes of volume L'^ one introduces the weight w^''{r) of the box 
of volume Lf^ around the point r 



wl'-ir) = [ d'^r'^l{r + r') 



(32) 
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that generalizes the microscopic weight ■0£ (r) corresponding to Lf, = 1 . Then all previous multifractal notions apply 
to the box weights if one replaces L by the ratio L/L^ of the two length scales. Equation [1] becomes 



Prob 



Lt 



and equivalently, their moments behave as 



da 



f{a)-d~qa 



L 



(33) 



(34) 



This formulation with coarse-grained boxes is also sometimes used numerically, in particular to measure correctly the 
scaling behaviors of negative moments q <0 [13, HI] . 



D. Analysis of the renormalization process 



To go from the microscopic scale / = 1 to the macroscopic scale Z = L of the whole system, we may introduce 
intermediate scales Im regularly placed on a logarithmic scale. For definiteness, let us consider a discrete system with 
L = 2^^ and introduce the intermediate scales 

/„ = 2" with m = 0, 1,..,A/ (35) 

We are interested into the flow of microscopic weight ilP'{r) upon renormalization. Up to now, we had always 
assumed the usual normalization of Eq. 2] on the full sample of volume L'^, but here to analyse the coarse-graining 
process using samples of various sizes, it is more convenient to work with fields '0^(r) free of any global normalization 
constraint. To characterize the microscopic weight ip'^{r) within the box of volume if-^ surrounding the point r, we 
introduce the variables 

If '0^(r + r') were constant within the box of volume if^, one would have W„i = 1. The lowest scale = 1 is 
characterized by Wq = ip'^ (r) / ip'^ (r) = 1. We are interested into the renormalization process 

Since the box of size l"^ contains the box of size ^m-ii it is convenient to introduce the positive exponent am > 

dVV-^lr + r') 

nam — '■m ^ (00\ 

~ ^ dVV2(r-Hr') ^ ^ 

to characterize the coarse-graining step from Im-i to Im- The ratio of two successive rescaled variables Wm of Eq. [351 
then reads 

w 2^L d^r'^P^r + r') 

in terms of the exponent 



The variable Wm can be then written as 



in terms of the accumulated value 



jm = am-d (40) 
Wm - 2-1^" (41) 



E^" (42) 
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along the RG trajectory. The multifractal definition for the weight 

Prob(-, K-T^-] « La^'''>d-f (43) 

in terms of the function .9(7) introduced in Eq. 1281 corresponds to the foUowing large deviation behavior for the 
variable rj\/ as M = In L/ In 2 becomes large 

Prob(—^j] cx 2^^»('') (44) 

Equivalently, the definition of the exponents A(g) from the weight moments at Anderson transition 
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c . ^ 1 « L-^(9' (45) 

corresponds, in the large deviation theory of the variable Tm, to the following statement for its generating function 

(2-«i^«) cx 2^^^'^(«) (46) 

The symmetry property of Eq. [29] for 5(7) means that the probability to obtain a value Tm ~ A/7 and the 
probability to obtain the opposite value Tm = —Mj are related via 

ProbjTM = iM) ^ 2A/[9(7)-9(-7)] = 2^'^ = 2^" f47) 
Prob{TM = -jM) ^ ' 

We now come to our central assumption. It seems physically natural to consider that the renormalization process 
of Eq. \3T\ represents some Markov chain, i.e. one expects that the probability to see Wm at scale 1^ will depend on 
Wni-i but will not depend on the previous values Wn with n < m — 1 that describe the finer structure at smaller 
scales inside the box of volume l^-i- The analogy with the analysis of non-equilibrium stochastic dynamics recalled 
in section IIII Al then suggests that the physical origin of Eq. [JT] is that the increment 7 of the considered functional 
Tm exactly measures the ratio of two opposite transition probabilities as in Eq. 1251 i-O- with our present notations 

r = (48) 

And indeed, using Eq. 1391 if oj^ic writes an elementary transition probability as 

k^{W -^W')=6{W' -2-'^W) (49) 

one obtains, using the properties of the ^-function, that Eq. [48] is satisfied 

k^{W W) = 2'<5{W - 2^M^') = 2T'fc_^(iy' ^ W) (50) 

This interpretation suggests that the symmetry relation of Eq. |^7| has an origin which is completely independent 
of the physical problem under consideration, so that it should be valid not only for Anderson transitions, but also 
for other random critical points. We should stress that for each random critical point, the large deviation function 
(7(7) is of course non-trivial and depends on the critical model, but the symmetry of Eq. [^5] seems generic and 
model- independent . 

The physical meaning of Eq. [JS] is that the variable 7 characterizes the irreversibility of the RG fiow. In a pure 
homogeneous system, the stability of the fixed point corresponds to 7 = 0, so that the forward RG and the backward 
RG are equivalent. In the presence of quenched disorder however, the exponent 7 becomes distributed, and its typical 
value jtyp corresponding to gijtyp) = is strictly positive 

Ityp > (51) 

This introduces some asymmetry between the forward RG flow and the backward RG flow. The strong multifractality 
limit where jtyp is far from zero corresponds to the far-from-equilibrium limit for dynamical systems, whereas the weak 
multifractality limit where ^typ is near zero corresponds to the close-to-cquilibrium limit for dynamical systems. Of 
course for a given random critical point, the disorder strength which has already been tuned to reach criticality cannot 
be tuned again to reach the regime of weak multifractality, in contrast with dynamical systems where one may always 
imagine to apply a weak external field to remain close-to-equilibrium. However, for a given type of random transition, 
there exists usually a parameter that allows to interpolate between weak multifractality and strong multifractality. 
For the usual Anderson transitions for instance, it is the space dimension d that leads to weak multifractality for 
d = 2 + e and to strong multifractality for large d [ll| . Similarly for two-dimensional random Potts models that will 
be considered in section [V] it is the number Q of possible states of an individual spin that allows to interpolate weak 
multifractality for Q = 2 -I- e and strong multifractality for large Q. 
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IV. SYMMETRY RELATION FOR MANY-BODY RANDOM CRITICAL POINTS 

111 the previous section, we have interpreted the symmetry of the multifractal spectrum discovered on the specific 
example of Anderson transitions as a Gallavotti-Cohen symmetry relation for a renormalization procedure. This 
suggests that this symmetry should be satisfied at other random critical points after an appropriate translation. In 
the remaining of this paper, we thus consider the case of many-body random critical points, like disordered spin 
models, where the multifractal properties are usually defined in terms of the statistics of the two-point correlation. 

A. Multifractal spectrum H{a) of the two-point correlation function 

Following the notations of Ludwig let us call (/)(r) the local order parameter, and X{N) the scaling dimensions 
of its disorder-averaged moments 0^ (?') . Then the two-point correlation function 

C{r,r + R) = {<P{r)(t){r + R)) (52) 

will present the following scaling behaviors for its disorder-averaged moments 

C^ir,r + R)^ J dec"" Pr{C) ^jx^j^^ (53) 

Equivalently, the probability distribution Pr{C) is described by the multifractal spectrum H{a) > such that 

dCPRic^-^] cx R-^^'-'^Ua (54) 

The saddle-point calculation of the integral in Eq. [53] yields that X{N) is the Legendre transform of H{a) 

- 2X{N) = max [-2Na - 2H{a)] (55) 

a 

The minimum cityp of H(a) where H[atyp) — governs the typical correlation 



In C(r, r + i?) = -2atyp In R (56) 
whereas the averaged correlation is governed by X{1) 
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B. Translation of the symmetry of f{a) into a symmetry relation for H{a) 

To translate the symmetry relation that holds for the singularity spectrum /(a) in Anderson localisation models to 
the many-body transitions, we should first rephrase the statements concerning two-point functions within an infinite 
system into statements concerning one-point functions within a finite-size system. In a system of size L'^, the local 
order parameter (t>L{i') is characterized by the scaling dimensions X{N) of Eq. [53] that govern the moments 

or by the multifractal spectrum H{a) of Eq. [54] 

d(t)PR ( (1) ^ ^] cx L-^'^'^Ua (59) 

\ L°- J L^ + oo 

Now to make the link with the Anderson localization framework or more generally with the usual multifractal for- 
malism, which is defined in terms of a normalized probability measure ll 
measure from the iion- normalized observables one is interested in [2^ . l44l . 
order parameter 



it is convenient to construct a probability 
45l . We thus introduce the normalized local 



(60) 
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where the denominator scahng is governed by the exponent X{1) of Eq. [58] as a consequence of the equivalence 
between spatial-averaging and disorder-averaging at the scaling level 

d'^rViW (61) 

The /(a) analogous to the Anderson localization definition of Eq. [T]is that the probability to have WL(r) ~ 
behaves as 



dwPR (wL{r) Lf^^'^-'^da 



(62) 



If WL{r) ^ the one-point function decays as (f>L{r) ~ 1/L°- with (see EgslSOland lCT]) 

a^a-d + X{l) (63) 

and thus the relation between the multifractal spectra /(a) and H{a) reads fEqslS^and lM)) 

f{a)-d=-H{a = a-d + X{l)) (64) 

Equivalently, the relation between their respective Legendre transforms r((j) and X{N) reads 

t{N) ^ X{N) ~ NX{l) + d{N ~l) (65) 

i.e. the exponents A(iV) introduced in Eq. [9] 

A(iV) = t(7V) - d(iV - 1) = X(iV) - iVX(l) (66) 

directly measure the non-linearity of X(N). 

The symmetry relation of Eq. [TUl translates for H{a) into 

H{2X{l)-a) = H{a) + a-X{l) (67) 

i.e. a symmetry with respect to — Xi. In terms of the Legendre transform X{N), the symmetry becomes 

X{N)-X{1-N) = {2N-1)X{1) (68) 

which is analogous to Eq. [TT] Finally, in terms of A(A^) of Eq. [SSI the symmetry is given by Eq. [T^j again 

A(iV) = A(l - iV) (69) 

V. TESTS OF THE SYMMETRY FOR THE 2D RANDOM POTTS MODEL 

The random many-body transition where multifractal properties have been the most studied, is the disordered 
two-dimensional ferromagnetic Potts model [l^, [l^, [ij, [l5|, [3, [13, [3 ■ it is a lattice spin model defined by the 
Hamiltonian 

-pH^Y. J^A„<r, (70) 

where the spin variables can take Q different values ct^ G {0, . . . ,Q — 1}. The sum extends over nearest neighbors on 
the lattice. The exchange couplings Jy are quenched variables. In the following, we have considered two different 
types of disorder realizations: i) self-dual random bond realizations, i.e. two different non zero 'ferromagnetic' values 
of the nearest neighbour interaction distributed with the same probability and ii) dilute bond, i.e. non zero values 
distributed over a fraction of the links only. The two types of disorder realizations belong to the same universality 
class. The phase transition is of second order for any value of Q in the disordered case (whereas in the pure case, the 
transition is of second order for Q < 4 and of first order for Q > 4). 

The moments of the spin-spin correlation function are expected to decay algebraically 



{S..^.r-\n~f,r^'''> (71) 

where we keep the prescription that . . . denotes the average over disorder while (. . .) is the thermal average. The 
exponents X{N) are the exponents that have been introduced in the previous section in Eq. [531 in the general case. 



10 



Numerical tests of the symmetry of the multifractal spectrum for various 2 < Q < 



We now present numerical data obtained using the connectivity transfer matrix introduced by Blote and Nightin- 
gale d^. In this formulation, even non- integer values of Q can be considered because the number Q of states enters 
only as a parameter. Furthermore large strips can be considered because the dimension of the transfer matrix is 
smaller in the connectivity space than in the spin space where it grows as x for strip width L. Finally the 
spin-spin correlation functions can be expressed in terms of connectivity and thus be computed very efficiently. In 
the following, the data have been averaged over 80,000 disorder realizations. Further details about the numerical 
procedure can be found in [l^ . 
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FIG. 1: Numerical test of the symmetry [69] for the quantity A{N), calculated from the exponents X{N) associated to the decay 
of the moments (ai<7i+„)^ for the self-dual random bond Potts model for Q < 3.5. The red symbols corresponds to the data 
after the operation N —f 1 — N while the black ones are the original data. The dashed line are the perturbative expansions: 
the first order of Eq. [77] is shown with short dashed lines, whereas the second order of Eq. [79] is shown long dashed lines. 



The system consists in an infinitely long strip of finite width L with periodic boundary conditions in the transverse 
direction, i.e. it is asymptotically one-dimensional. Therefore the correlations decay exponentially rather than alge- 
braically at the critical point and as a consequence of conformal invariance, the average correlation functions and its 
moments involve the critical exponents X{N) in a very simple manner: 



R) 



dCC^'PRiC) ^ e-T^W«, 



(72) 



In the strip geometry, the system thus exhibits an even closer relation to the dynamical model of section IIII Al with 
the role of time t now played by the coordinate along the strip. In figure [TJ we have plotted A(iV) obtained from the 
exponents X{N) using the definition Eq. [Mj The original data arc plotted as filled black symbols while the exponents 
A(A'^) after the transformation I ~ N —i- N are in open red symbols. Within numerical accuracy, the symbols fall on 
top of each other as expected from Eq. [§21 
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FIG. 2: Numerical test of the spectrum symmetry as calculated from the set of exponents X{N) obtained from the decay of 
the moments {(7icri+„)^ for the self-dual random bond Potts model for Q < 3. The filled circles correspond to H{a) vs a, while 
the open squares correspond to H{a) +a — X{1) vs 2X{1) — a. The part of the curves which superimpose (black and red online) 
is emphasized. 

The spectral function H(a) follows from the expression Eq. [72] of the spin-spin correlation function on the strip: 



H{a) = X{N) - aN. 



(73) 



Wc present on figure [2] the multifractal spectrum H{a) for the self-dual random-bond Potts model. This spectrum 
has been obtained using Eq. 1731 and the exponents X(N) computed by interpolation of the moments of the spin-spin 
correlation function with Eq. 1721 On top of the multifractal spectrum H{a), full symbols on the figure, we have 
plotted as open symbols the image of H{a) under the transformations H {a) + a — X{1) H{a) and 2X{1) — a ^ a. 
If the symmetry holds for the random-bond Potts model, the two curves should fall on top of each others according 
to Eq. 1671 We have emphasized with black and red colors on the figure the region where the two curves superimpose. 
The numerical data are in very good agreement with Eq. [67] for Q from 2.25 to 3. 

We now present another test of the symmetry using an alternative determination of the multifractal spectrum H{a). 
The definition 1721 written in terms of the variable y = — In C reads as a Laplace transform 



dYPR{Y)e 



-NY 



(74) 
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FIG. 3: Numerical test of the spectrum symmetry as calculated from the probability distribution of the correlation function 
of the dilute bond Potts model for Q > 3. The full lines (black online) correspond to H{a) vs a, while the filled circles (red 
online) correspond to H{a) + a — X{1) vs 2Jf (1) — a. 



The inverse Laplace transform evaluated close to the saddle point approximation leads to (see Ref. pTj ) 

P«(a)-(^)''%-^-(^) (75) 

with a = —IiiC/{2ttR/ L). The spectral fmiction is thus directly given in terms of the probability distribution 
Pr(— lnC/(27ri?/i)). Using again the connectivity transfer matrix, we have computed the probability distribution 
of the spin-spin correlation function for the diluted Potts model [131 ■ The multifractal spectrum was then obtained 
by interpolation of the numerical data with Eq. (75] In figure [31 the multifractal spectrum is plotted in black 
for three different values of the number of states Q. In red, we have plotted the data after the transformations 
H {a) + a — X{1) H{a) and 2X(1) — a ^ a. Again, the data fall nicely on top of each other for small and 
intermediate values of a. For large values of a, the two curves significantly deviate from each other, but at the 
same time, the data dispersion which measures the uncertainty also increases. According to 1731 this region, where the 
multifractal spectrum increases with a, is associated to negative moments N of the spin-spin correlation function. The 
latter are dominated by rare events corresponding to weak spin-spin correlation and thus strongly diluted realisations 
of disorder. Numerically, one can hope to sample correctly events whose probability is of order 1/A^Dis. where Nuis. is 
the number of disorder realizations considered. Very rare events arc unlikely to be generated. Wc may thus interpret 
the deviation observed on figurc[3las due to an underestimate of the tail of the probabihty distribution P{C) < 10^^. 
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B. Reminder on existing perturbative expansions in (Q — 2) 

The Harris criterion indicates that disorder is relevant at the transition for Q > 2 but only marginal for Q = 2: the 
point Q = 2 can be thus used as a starting point for a perturbative expansion in the parameter (which measures the 
relevance of disorder) 

VH = a/z. = i-(Q - 2) - ^(Q - 2)2 + 0{Q - 2f (76) 
which measures the disorder relevance [131. 



1. First order contribution 

The exponents X{N) governing the correlation functions have been computed at first order in yn by Ludwig [l3 |: 

XiN) = NXil) - ^^^^^^\ h + 0{yl) (77) 

The corresponding Gaussian multifractal spectrum H{a) can be foimd in Ludwig [l3 |. but the symmetry is much 
easily tested under the form of Eq. [53 According to the definition Eq. [511 the expansion Eq. [77] leads to 

A(iV) = XiN) - NX{1) ^ - ~ + Oiyj,) (78) 
which is manifestly invariant under the transformation N ^ I — N . 



2. Second order contribution 



The perturbative expansion of Ludwig has been extended to the second order in a replica symmetric scenario by 
Lewis p^ : 



X{N) = NX{1) 



N{N -1) 
16 



11 , N -2 f 29n 
IJH + I — - 41n2 + — ( 33 - 
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24 



V3 



(79) 



The presence of the (TV — 2) term explicitly breaks the symmetry iV ^ 1 — iV. As a comparison, the equivalent 
weak-multifractality perturbative expansion for Anderson localization can be found in [29j up to order in e = d — 2 
and satisfy the multifractal symmetry. 

If the result of Eq. [79|is really true, this means that the multifractal symmetry is not realized in the random Potts 
model. However, we believe that one should be cautious concerning the result of Eq. 1791 in spite of previous numerical 
confirmations. It is a perturbative replica result that has been derived within the so-called replica-symmetric scenario. 
Indeed in some disordered systems, the replica symmetry is known to be spontaneously broken, so that one should in 
principle study the various fixed points corresponding to various replica symmetry properties and then decide which 
one corresponds to the physical fixed point. For the random Potts model, two different types of fixed points have 
been compared in [ZJ: the replica symmetric fixed point, used to derive Eq. 1791 and the broken replica symmetry 
fixed point based on the Parisi block diagonal matrix Ansatz. We refer to [i^ for a detailed discussion of the possible 
physical meanings of replica symmetry breaking at the critical point of random ferromagnets. According to [47|, the 
differences in the critical exponents between these two fixed points appear only at second order y\^ so that the first 
order contribution of Eq. [771 although derived within the replica symmetric fixed point, seems to be of larger validity 
in the space of replica possible fixed points, whereas the second order contribution of Eq. [79l seems specific to the 
replica symmetric fixed point. Since it has not been possible up to now to determine on theoretical grounds which 
fixed point describes the scaling limit of the disordered lattice models of Eq. I70[ very precise numerical computations 
16 , 4£ I have been necessary to discriminate between the two fixed points studied in [47| : these numerical studies 
16l . |49| have concluded (i) that the replica symmetric broken fixed point based on the Parisi block diagonal matrix 
Ansatz is excluded, and (ii) that the replica symmetric fixed point is thus the physical one. However we believe 
that (i) does not directly imply (ii), because only two particular fixed points have been studied. Since the Parisi 
block diagonal matrix Ansatz comes from spin-glasses in the mean-field limit, whose physics is completely different 
from two-dimensional random ferromagnets, one may argue that another type of replica symmetry breaking could be 
realized in critical random ferromagnets. For instance, it has been proposed in [sot that in some disordered systems, it 
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was appropriate to consider some vector symmetry breaking scheme instead of the usual Parisi matrix sector breaking 
Ansatz. Very recently, still another pattern of replica symmetry breaking has been identified to account for a pre- 
freezing phenomenon of multifractal exponents in a disordered system [sij. As a consequence, the proper identification 
of the true physical fixed point among the possible fixed points in replica space remains the important issue to decide 
whether Eq. [79] is really true and sufficient to rule out the multifractal symmetry. 

C. Discussion 

In summary, we cannot give a definite conclusion on the existence of the multifractal symmetry in the random Potts 
model, but we would like to stress the two following points: 

(1) from a purely numerical point of view, the data presented on the Figures above satisfy the symmetry within 
the error bars for all values 2 < Q < % that have been studied. This statement is valid for the self dual disordered 
model and in the dilute case. Moreover, the differences between the numerical data and the second-order perturbative 
expansion of Eq. [75] plotted as dashed lines in figure [I] are significantly larger than the deviation from a symmetric 
A(iV). 

(2) from a theoretical point of view, it would be interesting to study other fixed points that break the replica 
symmetry using various schemes, to see whether one can identify a type of scheme that would preserve the multifractal 
symmetry. 



VI. CONCLUSIONS AND PERSPECTIVES 

In this paper, we have analyzed the physical origin of the symmetry relation f{2d — a) = f{a) + d — a proposed by 
Mirlin, Fyodorov, Mildenberger and Evers [27| for the singularity spectrum /(a) of critical eigenfunctions at Anderson 
transitions. We have explained the analogy with the Gallavotti-Cohen symmetry of large deviation functions that are 
well-known in the field of non-equilibrium dynamics: the multifractal spectrum of the disordered model corresponds 
to the large deviation function of the rescaling exponent 7 = {a — d) for a renormalization procedure considered 
as a Markov chain in the 'effective time' t = InL. We have concluded that the symmetry discovered in the specific 
example of Anderson transitions should actually be satisfied at many other random critical points after an appropriate 
translation. For many-body random phase transitions, where the critical properties are usually analyzed in terms 
of the multifractal spectrum H{a) and of the moments exponents X(N) of two-point correlation function, we have 
obtained that the symmetry becomes H{2Xi — a) = H{a) + a ~ X{1), or equivalently A(iV) — A(l — N) for the 
anomalous parts A(iV) = X{N) — NX{1). We have presented detailed numerical tests in favor of these relations, 
or at least compatible with them, for the two-dimensional random Potts model with various Q > 2. Although 
presently available numerical results for the scaling dimensions are in favour of the RG results deduced from the 
replica symmetric scenario (which does not satisfy the multifractal symmetry) and rule out the Parisi matrix sector 
breaking Ansatz, we cannot exclude that other fixed points would equally be compatible with the numerical results, 
but also with an exact multifractal symmetry. It would be interesting in the future to test this multifractal symmetry 
at other random critical points. At a more general level, we hope that the analogy with the field of non-equilibrium 
dynamics will give new insights into the properties of renormalization fiows in disordered systems. 
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